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The features of equatorial motion of an extended body in Kerr spacetime are investigated in 
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I. INTRODUCTION 

The description of the gravitational interaction between the constituents of a binary system in the general theory of 
relativity requires taking into due account their internal structure. The orbital dynamics of two bound compact objects 
is tackled in the literature by a plenty of different methods resorting to various approximation schemes. Analytic 
approaches include notably the post-Newtonian approximation [1], possibly implemented using effective field theory 
techniques [2], as well as the gravitational self-force corrections to geodesic motion [3], which can both be combined 
efficiently by means of the “effective-one-body” approach [4, 5]. 

When the mass of one body is much smaller than the other, the problem boils down to studying the dynamics 
of an extended body in a fixed background field, generated by the heavier mass. In this approximation, a self- 
consistent model describing the evolution of both linear and angular momenta for pole-dipole sources was developed 
by Mathisson [6], Papapetrou [7, 8], Pirani [9], Tulczyjew [10], and later generalized to bodies endowed with higher 
multipoles by Dixon [11-15]. The Mathisson-Papapetrou-Dixon (MPD) model accounts for the motion, on a fixed 
background, of a point-size test object with internal degrees of freedom, in the absence of significant gravitational 
back reaction. 

The main astrophysical situation for which a full relativistic treatment is needed occurs when the object — assumed 
to be compact to avoid tidal disruption — experiences the strong field produced by a nearby black hole. In that case, 
the MPD approach may be used to investigate the evolution of the system, but the parameters of the model must 
be regarded as effective ones [16]. Although self-force effects are not negligible on larger-than-orbital time scales, 
they only yield higher-order corrections. On the other hand, combined with dimensional regularization, the MPD 
model is appropriate to describe the dynamics of (self-gravitating) compact binaries including finite-size effects, 
in the post-Newtonian framework [17-19]. It yields results that are dynamically equivalent to those derived from 
suitable effective actions [20-22]. The “skeleton” stress-energy tensor encoding the MPD evolution is known at the 
quadrupolar level [23]. Its octupolar contributions have also been obtained recently [24] assuming an effective, Bailey 
& Israel type, Lagrangian [25] . Those corresponding explicit expressions have been used to build accurate theoretical 
templates for the signal of gravitational radiation emitted by those sources [19, 24, 26, 27], in the context of the 
data analysis of gravitational-wave observatories, such as the advanced Virgo [28] and LIGO [29] detectors, the future 
cryogenic interferometer KAGRA [30] or, possibly, the space-based observatory eLISA [31], a candidate for the future 
L3 mission of the European Space Agency. Most of post-Newtonian expressions can be checked by comparing them 
to the test-body counterparts, in the extreme mass-ratio limit. 

In this paper we study the dynamics of an extended body endowed with both spin and quadrupole moment in a Kerr 
spacetime using the MPD model. The motion is assumed to be confined on the equatorial plane, the spin vector of 
the body being aligned with the axis of rotation of the central object. In previous works, we have discussed the effects 
on the dynamics of a general quadrupole tensor in both Schwarzschild and Kerr spacetimes [32-34]. Here, we consider 
more specifically the case of a spin-induced quadrupole tensor. We assume that the object reaches thermodynamic 
equilibrium in its proper frame on time scales that are very short compared with the orbital period and neglect the 
tidal deformations. Its internal state thus depends adiabatically on the mass and the (instantaneous) spin. Using 
effective field theory arguments, it is then straightforward to check that the body quadrupole is actually quadratic in 
the spin. This situation was described in details by Steinhoff & Puetzfeld [35], who developed a very general framework 
to include quadratic in spin corrections as well as tidal interactions in the MPD scheme, with special attention to the 
study of the binding energy of the system as obtained from the analysis of the associated effective potential. Later, 
Hinderer et al. [36] performed an analysis of the corresponding dynamics, in order to compare the periastron advance 
and precession frequencies with those of a different approach, but restricted themselves to a very special (although 
physically motivated) choice of the quadrupole tensor, leading to great simplifications in the analytic treatment. 

Here, we shall assume the same form of the quadrupole tensor as in Ref. [35], but neglect quadrupolar tides, i.e., 
our quadrupole tensor is of the electric-type only and is proportional to the trace-free part of the square of the spin 
tensor by a constant parameter, whose numerical value is a property of the body under consideration. For neutron 
stars such a quantity depends on the equation of state [37], while it is exactly 1 for black holes. We keep it a free 
parameter of the model that can affect associated observables, like the energy and the angular momentum, which we 
computed explicitly and compared with the results of Ref. [36], or the Innermost Stable Gircular Orbit (ISCO) and its 
frequency, discussed here in details. We achieve a fully analytic treatment of the MPD equations in a “perturbative” 
scheme, obtaining corrections to geodesic motion up to the second order in spin. 

Throughout this work we use geometrical units, setting the Newton constant G and the speed of light c to 1. Tensors 
are represented either in abstract notation or in index notation combined with the Einstein convention, depending 
on the context. Greek indices refer to spacetime coordinates and vary from 0 to 3, i.e., /r = 0,1, 2,3, whereas Latin 
indices, ranging from 1 to 3, label space coordinates. The notation stands indifferently for the partial derivative 
with respect to the ^th coordinate or for the coordinate basis vector associated to x^, while d denotes the exterior 
derivative. The spacetime metric g^i,, taken to have signature defines a unique Levi-Givita covariant 
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derivative, and an associated Riemann curvature with the convention that = (VqV /3 — V/ 3 Vq,)u^ 

for any vector field . Symmetrization of a tensor T over a set indices is indicated by round brackets enclosing them: 

_j_ T''>^)/2. Instead, for index antisymmetrization, square brackets are used: — T^^‘)/2. 


II. MPD DESCRIPTION OF EXTENDED BODIES 


In the quadrupole approximation, the MPD equations read 


dpm 

dr 


= -L 

(spin) ~ (quad) ’ 


1 


dr 


2plt^U''^ + 1 
3 


(spin) (quad) 


( 2 . 1 ) 


( 2 . 2 ) 


where = mu^ (with u - u = = —1) is the total 4-momentum of the body with mass m and direction u^, 

is the (antisymmetric) spin tensor, is the quadrupole tensor, and = dz^/dr is the timelike unit tangent 

vector — or 4-velocity — of the body “reference” line, parametrized by the proper time r (with parametric equations 
= z^(t)), used to make the multipole reduction. 

In order to ensure that the model is mathematically self-consistent, the reference point in the object should be 
specified by imposing some additional conditions. Here we shall take the Tulczyjew conditions [10, 11], 


= 0. 


(2.3) 


With this choice, the spin tensor can be fully represented by a spatial vector (with respect to u), 

s(ur = irfiur^^s^'y = , ( 2 . 4 ) 

where r](u)a/ 3 'y = is the spatial unit volume 3-form (with respect to u) built from the 4-volume form rja/B'YS = 

V~9 £c(/ 37 ( 5 , with Cap-yS (eoi 23 = 1) being the Levi-Civita alternating symbol and g the determinant of the metric in 
a generic coordinate grid. Using a fairly standard convention, hereafter we denote the spacetime dual of a tensor S 
(such that /2) by *S, whereas the spatial dual of a spatial tensor S with respect to u is represented 

by *(“>5'. It is also useful to introduce the signed magnitude s of the spin vector, which is not constant in general 
along the trajectory of the extended body: 

= S{ufS{u)f> = = -iTr[52], (2.5) 

with Tr[S'2] = 


A. Spin-induced quadrupole tensor 

The 1-1-3 decomposition of the quadrupole tensor J and its general properties are briefly reviewed in Appendix A. 
We will consider here the physically relevant case where it is completely determined by the instantaneous spin structure 
of the body (see, e.g.. Refs. [35, 38]). More specifically, we shall let the quadrupole tensor have the form 

ja/375 = 4u[“A(u)^]['>'lt^] , (2.6) 

with 

= (2,7) 

4 m 

where Cq is a “polarizability” constant and [5'^]®"'"^ denotes the trace-free part of , i.e., in terms of both the 

spin vector and the associated spin invariant, 

[S'2]STFa/3 ^ _ \p{u)°‘PSp^S'^P = S{u)°^S{uf - \s‘^P{uT^ . 

o o 


( 2 . 8 ) 
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The values of Cq associated with compact objects are given, e.g., in Ref. [39]. The normalization is such that Cq = 1 
in the case of a black hole [40], whereas for neutron stars Cq depends on the equation of state and varies roughly 
between 4 and 8 [37]. 

For the spin-induced quadrupole tensor (2.6), the link between and takes a particularly simple form if cubic- 
in-spin corrections are neglected. Indeed, contraction of Eq. (2.2) with Pn shows that, apart from corrections of order 
J = 0{S^), the difference (P • U)P^ — P'^U^ is precisely Pi,DS^’^/dr = —S^^DP^/dr, which can be seen to be of 
order 0{S^) from Eq. (2.1). In the end, P^ is approximately proportional to and, as an important consequence, 
the right-hand side of the precession equation (2.2) is at least quadratic in the spin. As, on the other hand, the 
right-hand side of the precession equation (2.2) is at least linear in the spin, we conclude that the time differentiation 
of our kinematical variables actually multiply any combination of them by a factor 0{S). 

In Ref. [36], the MPD description of test bodies endowed with a spin-induced quadrupolar structure is used to 
check the consistency of the computation of the periastron advance for a binary system with the effective-one-body 
formalism in the extreme mass-ratio limit, including terms that are quadratic in the spin. The quadrupole tensor is 
assumed there to be that of a black hole and to take the form (2.6) with Cq = 1, so that we must recover the results 
of Ref. [36] for this particular value. 


B. Simplified form of the MPD equations 

The MPD equations can be written in a more convenient form at quadratic order in the spin. First, the quadrupolar 
contribution in Eq. (2.1) splits into a parallel and a perpendicular part with respect to the direction 

^r^uad) = [P{U)r^.Ro.MS + ^ ■ (2.9) 

Next, neglecting remainders that are cubic in the spin, the quadrupole and the velocity in the last term may 

be moved under the operator D/dr, since their covariant time differentiation would actually produce terms smaller 
than the original ones by a factor 0{S), as explained in the previous subsection. The equations of motion then 
become: 

D 1 1 D 

^ = - 2 U" - g [P{U)rV,R^p^s + ^ {mjU^^) , (2.10) 


where we have posed 


( 2 . 11 ) 

Thus, defining a modified linear momentum = P^ — mjU^ effectively changes the quadrupolar force into 

its projection orthogonal to C/^, while the spin force is orthogonal to U^. The precession 

equations are unaffected. More explicitly, Eqs. (2.1) and (2.2) expressed in terms of become: 

Dn^ 1 1 

^ = -r 1/" V^Ro^p^s + 0(^3), (2.12a) 

dr z D 

D 4 

—— = 2p[^C/'^l -f - -b 0(S'3) , 

dr 3 

with = [P{U)Y'^S/v Now, if we multiply the equations of motion (2.12a) by Pp oc t/p -I- 0{S^), we see that, at 
our accuracy level, the first term on the right-hand side vanishes by virtue of the Riemann-tensor symmetries whereas 
the second one is zero due to the contraction of t/^ with [P(t/)]^'". We conclude that the effective mass 
is conserved, modulo cubic spin corrections (see also Appendix B). It may be regarded as the “bare” mass of the 
extended body, mg, at the quadratic order in the spin, so that its mass m to order 0(5'^) is given by 

m = mg + mj + O(S^). (2.13) 

Finally, the MPD equations (2.1)-(2.3) imply that the unit vectors U and u are related by 

+ 0{S^) . 


(2.14) 
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C. Conserved quantities 

In stationary and axisymmetric spacetimes endowed with Killing symmetries, the energy E and the total angular 
momentum J are conserved quantities along the motion, associated with the timelike Killing vector ^ = dt and the 
azimuthal Killing vector rj = respectively. They are given by 

E = -e„P“ + , 

J = , (2.15) 

where 

pfp = = 94>[a,P] , (2.16) 

are the Papapetrou fields associated with the Killing vectors. Note that E and J as defined above are conserved to 
all multipole orders in spite of the higher multipolar structure of the body, which is entirely encoded in P [41]. 

The conserved quantities (2.15) for a purely dipolar particle in a Kerr spacetime have been computed, e.g., in 
Refs. [42, 43]. The expressions given there are general enough to account for all higher-order spin corrections but 
those coming from the spin-induced multipole moments (i.e., are exact when the quadrupole and higher multipole 
moments vanish). This means in practice that our results should reduce to those of Refs. [12, 43] for Cq = 0. 


III. MOTION IN A KERR SPACETIME 


The Kerr metric in standard Boyer-Lindquist coordinates (t, r, 0, </>) reads 
2Mr\ 


ds^ = - 1 - 


, 2 4aMr . o ^ i 2 ^ (?'^ + a^)^ — Aa^ sin^ . 2 1 ,2 1 \ 

dt^-——sin^ 0dt d^-I-— dr^-I-Sd0^-I--—-sin^0d0 , (3.1) 

2-/ ZA 2-i 


with A = — 2Mr + a^ and E = r'^+a^ cos^ 6. Here, a > 0 and M > a denote the specific angular momentum and the 

total mass of the spacetime solution, respectively. The event and inner horizons are located at = M ± y /— a^. 
Let us introduce the zero angular-momentum observer (ZAMO) family of fiducial observers, with 4-velocity 

n = N-\dt-N‘^d^) (3.2) 


orthogonal to the hypersurfaces of constant t, where N = (— 3 “) and = gt,p/gci,ci, are the lapse and shift 
functions, respectively. A suitable spatial orthonormal frame adapted to ZAMOs is given by 


with dual 


— Tl ^ — 


\/ 9rr 


dr — d-r , Gq — 


— da 1 G2 — _ d(p — dl , 




■'e ’ ^((> 




u}^=Ndt, U}'^ = y/^dr, uj^ = ,/^d6, = ^/^{d(f> + N^^dt). 


(3.3) 


(3.4) 


The ZAMOs are subject to the acceleration a{n) = V„n. They are locally non-rotating, in the sense that their 
vorticity vector uj{n)°‘ vanishes due to their surface-orthogonal character, but they have a nonzero trace-free expansion 
tensor 9[n)ap = P{n)^^P(nY the latter, in turn, can be completely described by an expansion vector 
O^inY = 0 (n )“/3 such that 


e{n) =e^® e^{n) + e^{n) ® 


(3.5) 


where ® represents the tensor product. The nonzero ZAMO kinematical quantities (i.e., acceleration and expansion) 
all belong to the r-0 2-plane of the tangent space [44-47], with 

a(n) = a{nYCf + a{nYeg = df{hiN)ef -f clg(ln7V)eg , 

(9,A^e, + dgN^eg) . 


2N 


(3.6) 
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It is also useful to introduce the curvature vectors associated with the diagonal metric coefficients, 

k{x'‘ ,n) = k{x^ ^nYef + n)®eg 

= - [dfYii y/g~)ef + dY^iiy/g~)eY\ ■ 


(3.7) 


We shall use the notation K{cj),ny = fc(Lie) for the Lie relative curvature [45, 46], largely adopted in the literature, 
and limit our analysis to the equatorial plane {0 = tt/2) of the Kerr solution, where 


N = 


rA 


1/2 


+ a^r + 2a?M 


N‘l> 


2aM 

r3 _|_ fj2j, _|_ 20^M ’ 


(3.8) 


and A = The ZAMO kinematical quantities as well as the nonvanishing frame components of the Riemann 

tensor are listed in Appendix C. 

Let us now consider a test body rotating in the equatorial plane around the central source. Its 4-velocity U may be 
written in terms of the velocity n) = v'^ef + relative to the ZAMOs, with associated Lorentz factor 7 (C/, n), 
as 

[/= 7 (C/,n)[n + z/(C/,n)], 7 ( 1 /,n) = (l - ||i/([/,n)|p) , (3.9) 

The parametric equations of the orbit are solutions of the evolution equations U = da;“/dr, i.e.. 


dt 7 dr 71 /’" 

dr A ’ dr y/grr ’ 

^ ^ 

dr V N J ’ 


(3.10) 


where the abbreviated notation 7 ( 1 /, n) = 7 and = iy(U,nY has been used. For equatorial orbits, a convenient 
parametrization can be r itself instead of the proper time r. 

A case of particular importance is that of uniform, circular equatorial motion. The unit tangent vector, U, may then 
be parametrized either by the (constant) angular velocity with respect to infinity C or equivalently by the (constant) 
linear velocity v with respect to the ZAMOs, i.e., 

t/ = r[5t + C50] =7[e£ + z/e^], y = (1 - i/2)-i/2 ^ ^3^^^ 

with 

r = [N^ - g^^iC + = ^ , C = -A^ + ^^. (3.12) 

The parametric equations of the orbit reduce to 

t = to+TT, r = ro, ^ = (j) = 4>o + ilr, (3.13) 

where O = T^ is the proper time orbital angular velocity. 

For timelike circular geodesics on the equatorial plane, the expressions of the angular and linear velocities do depend 
on whether the orbits are co-rotating (+) or counter-rotating (—). They read 


respectively, with = \/M/r^ denoting the Keplerian angular velocity for a non-spinning, Schwarzschild black hole. 
The Lorentz factor of the corresponding 4-velocity U± is found to be 


r± 


IC±I 



3M 


± 2 aCif 


-1/2 


(3.15) 


In the static case, we can actually use the Schwarzschild values C± —>■ iC/f a-nd i'± —>■ with = y/M/{r — 2M). 

It is convenient [see Eq. (4.21)] to introduce a spacelike unit vector U± that is orthogonal to U± within the Killing 
2 -plane, by defining 

= f±[(9t -k C±9,/] = ±j±[iy±ef + e^] , 


(3.16) 
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in terms of the parameters 

f± = r±|:/±| = |ri±|-^ + — ip 2aCK^ , 

r I ^ 1 _ vCk I-2M/r±aCK 

^ y/g^ M 1 + a2/r2 =F 2aC_ft: ’ 


(3.17) 


where fl± = r±C± and where the ± signs correlate with those of U±. Note that C± = E±/L± is the ratio between 
the energy E± and the azimuthal angular momentum L± per unit mass of the particle. Those two quantities are 
expressed as 


E± 


L± 


N^± 1 + 


2aM 


r-\/A 


= 


1 ^ 

Ck 


l- — ±aCK 
r 


1 + ^ T 2aCif 


(3.18) 


A. Orbit of the extended body 

In order to describe the motion of the extended body according to the MPD model, we need both the timelike 
unit vector U tangent to the center world line and the unit timelike vector u aligned with the 4-momentum. In the 
following, we shall assume that the world line of the extended body is confined onto the equatorial plane, so that the 
4-velocity U is given by Eq. (3.9), with 

v{U,n) = vv = tf -I- = i^(cosaef -I- sinae^), (3.19) 

where v and a. € [0, are the signed magnitude of the spatial velocity and its polar angle, measured clockwise from 
the positive (j) direction in the r-(j) tangent plane, respectively, while u = v{U, n) is the associated unit vector; hence, v 
has positive/negative values for co/counter-rotating azimuthal motion (a = 7r/2) and outward/inward radial motion 
{a = 0) with respect to the ZAMOs, respectively. 

A similar decomposition holds for the (body) 4-momentum P = mu, in the case of equatorial orbits: 

u = 7 „[n -I- VuVu ], 7 „ = (1 - , (3.20) 

with 

v(u, n) = Vu = cos ttuCf -f sin , (3.21) 

and au G [0, ^]. An orthonormal frame adapted to u = Cq can then be built by introducing the spatial triad: 

ei = = sina„ef - cosQ;„e^ , 62 = sgn(t'„) 7 u(i/un-|-z)^), 63 =-eg. (3.22) 

The dual frame of {ca} will be referred to as {w“}, with w® = —u^, being the covariant dual of u. The projection 
of the spin vector into the local rest space of u defines the spin vector S{u) (hereafter simply denoted by S, for short). 
In the frame (3.22), the spin S decomposes as 

S = S^ei + 8^62 + 8 ^ 63 . (3.23) 


B. Setting the body’s spin and quadrupole in the aligned case 

In the following, we shall consider the special case where the spin vector is aligned with the spacetime rotation axis, 
i.e., 

8 = 863 . (3.24) 

This entails that the spin and quadrupole terms entering the right-hand sides of Eqs. (2.1) and (2.2) decompose, with 
respect to the frame adapted to m, as 


^(spin) — -P’(spm)^ + -^(Ipin)^! + ^(spin)e2 






and 


with 


^(quad) — -f"(quad)^ + -^(quad)®l + -^(quad )®2 , 

-^(spin) ^ ^ ^(spin)(^) 5 

-^(quad) ~ ^ ^ ^(quad) ('^) ? 


^(spin) ('^) — ^(spin)l^ “t“ ^(spin) 2 d^ : 

^(quad)(^) ~ ^(quad)l^ “t” ^(quad)2d^ 5 

respectively. The explicit expressions for the above components are listed in Appendix D. 


(3.25) 


(3.26) 


(3.27) 


IV. SOLVING THE MPD EQUATIONS FOR NON-PRECESSING EQUATORIAL ORBITS 

A. Complete set of evolution equations 

Under the assumptions of equatorial motion and aligned spins discussed in the previous section the whole set of 
MPD equations (2.1)“(2.3) reduces to 

^ _ po , pO 

'^(spin) + ^(quad) 1 


=- — [pcos(a„ +a) - Vu] + — (sin aua(n)’' + sina fc(Lie)) 

1 


m-iuVu 
di^u 7 


(-^(spin) + '^(quad) ) ; 


dr 

ds 

dr 


Iti 


cos auttiriY + V sin(au + a)0i (n)' 




j(^(spin) + ^(quad)) > 


= 0 , 


together with the compatibility conditions 


0 63 ^ (^(spin) (^) “f ^(quad) (^)) “f ^ (-^(spin) “f -^(quad)) 


(4.1) 


(4.2) 


which come from the spin evolution equations and yield two algebraic relations for v and a. The last equation (4.1) 
implies that the signed spin magnitude s is a constant of motion. Finally, Eqs. (4.1) must be coupled with the 
decomposition (3.9) and (3.19) oiU = da;“/dT to provide the remaining unknowns t, r and </> [see also Eqs. (3.10)]. 
Eqs. (4.1) and (3.10) are written in a form that is suitable for the numerical integration. Additional (non-independent) 
relations are obtained from the conservation (2.15) of the total energy and angular momentum, 


E = Nju 


m + s{vu sinaua{nY + S^{nY) 


-N^J, 


J = luy/^ miyuSinau - s(fc(Lie) + sina„6>^(n)’') 


(4.3) 


These can be used as a consistency check. Examples of numerically-integrated orbits are discussed in Refs. [33, 34]. 

We are interested here in studying the general features of equatorial motion to second order in spin, taking advantage 
of the simplified form of the MPD equations discussed in Section IIB. We shall also derive analytic solutions for the 
orbits to that order by computing the corrections, produced by a non-zero spin, to a reference circular geodesic motion. 


B. Solution to the order O(S^) 

Let us introduce the dimensionless spin parameter 


s = 


moM 


(4.4) 
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where mo denotes the “bare” mass of the body. We shall systematically neglect terms that are of orders higher than 
the second in s, hereafter. Hence, all quantities must be understood as being evaluated up to the order 0{s^). Our 
set of equations can then be simplified by means of Eqs. (2.13) and (2.14), which yield 

m = mo - icgmoM^s^y^ [{2Eff + Egg)vl cos2a„ + (2 + vl)Egg + Avu sinauH-g] , (4.5) 

and 

V = Vu + ]^{^- Cq) [vu [(2Eff + Egg) cos2au + 2>Egg] + 2H~g{l + vl) sina^} , 
a = a„ - (1 - Cq) M'^s^ [{2Eff + Egg)vu sina„ - H.g] , (4.6) 


respectively. Here, Eff, Egg are components of the electric part of the Riemann tensor, while H~g is a component 
of its magnetic part (see Appendix C). Eqs. (4.6) show that the value of the polarizability for black holes, namely 
Cq = 1, is a very special case leading to a great simplification (see below). 

On the other hand, solving Eqs. (4.3) algebraically for and au leads to 


= 


7„j^„sina„ = 


E + N'^J 
rngN 


1-Ms 
1 + Ms 


hinr 


Na{nYJ 


{E + mj)^ 


^"^'(moM2 + ’ 


e.{nf+ {E + N^J) 


y/ 94>4>^0-‘ie) 


NJ 


- 


f m2 


\mQM^ 


“t” E^'P “t” 


(4.7) 


Those identities can be used next at the lowest order in the previous equations so as to express m, u and a in terms 
of E and J. In particular, the solution for the body mass is found to be 


m = mo<jl--CQ^ 


1 + 


2 9 

rrinV^ 


{J-aEf 


s^ > = mo + s^m2 . 


(4.8) 


Its behavior as a function of the radial coordinate is shown in Fig. 1 for selected values of the parameters. It is 
interesting to evaluate the difference between the limiting values at the horizon and at infinity, |m(+oo) — m(r_|_)| = 
s 2 m 2 (r+), since this represents the largest mass variation during the evolution. 

It can be verified that the solution (4.6) for U is equivalently obtained from the compatibility conditions (4.2). 
Furthermore, the truncated equation (4.1) for the mass reads 


with 


^(quad) 


+ 0(s3), 


— 

(quad) 


-moM'^s^CQ-fuVu cos 


Wu{bi cos 2 a„ + C 3 ) + 2 (ai - 64 )sina„] + i( 2 ci 



(4.9) 


(4.10) 


where the lowest order piece of the solution (4.7) may be used, which implies notably that ^"(spin) = C(s^). The 
quantities Oi, and Ci, all functions of r, are listed in Appendix D. Once parametrized with the radial coordinate r 
instead of the proper time r by means of Eqs. (3.10), formula (4.9) takes the very simple form 


dm 2 

dr 
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(J-aE)^ , 


(4.11) 


whose solution coincides with that of Eq. (4.8). 


C. Circular motion 

i. Solution to the order 0{s^) 

In this subsection, we restrict ourselves to circular orbits, as described by the parametrizations (3.11)-(3.13). For 
circular motion in the equatorial plane, we must set a = 12, so that Eqs. (4.6) become 


V — Vu E — Cq) M'^S^ \—V^(^Ef.f. — Egg) + H~g{l + iz^)] , 
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r 


Figure 1. The behavior of the mass m of the body (in units of mo) as a function of the radial coordinate (in units of M) is 
shown for the following choice of parameters: a/M = 0.5, Cq = 1, s = 0.25, E/mo = 1 and J/{moM) = 4. The mass shift 
in this case equates |m(+oo) — m(r+)| « 0.055mo. The value of the dimensionless spin parameter s has been exaggerated in 
order to enhance the effect. 


TT 

““ “ 2 ’ 


(4.12) 


to second order in s. Thus, ^ from its expression (4.10), and Eq. (4.9) tells us that the mass of the body is a 

constant of motion to that order. However, it differs from the bare mass toq by virtue of the general definition (4.5), 
which yields 


m = mo - -CgmoM^s^y^ [Egg - + 2v±H-g] 


= mo 


^ 1 ^ =F Aay/Mr + 3a^ ^2 

1 rV 2 ^3/2 _ 3^^ ± 2av^^ 


once evaluated at = 7r/2, with replaced by iy± to the lowest order. 

Next, solving algebraically the second equation (4.1) for we find 

l^u = + Sl^u(l) + S^^u{2) 1 

where v± denote the geodesic linear velocities (3.14), while z^„(i) and z.^„( 2 ) are the spin-induced corrections 


M 

^u(l) = i [~^±{Efr — Egg) -|- (1 
^u{2) — (^^Mquad) 


'±)Hf.g] , 


'Ml) 


'■'u{l) 




«(Lie) 


1 


2Cic 

of first and second orders, respectively, with (see also Eq. (E5)) 


H-± (1 - Cq) MC,k ± 


MH. 


fe 


‘2'1±v±C.k_ 


F, 


1 


(quad) 


— 7^Q \ 7 ±[^i ~ ^2 + — ^5)] — + 7(^2 + ^3) f ■ 


Substituting the above solutions for and into Eq. (4.3) we obtain 

E = moNj± I 1- — - 12± I 

7±C±^^„(i)V^ + MiV(z.±a(n)^ + 0Hn)O + M^N^{v±eAny + 


mo7± 


(4.13) 


(4.14) 


(4.15) 


(4.16) 
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N 


■ mo7± 


^/9^C± + '^u(2) + 27±^±^^u(1) 




J = nno')±v±^/^ 

+ TOo 7 ±V^ 1 ±^U( 1 ) - M{v±e^{ny + A:(Lie)) 


m2'y±v±y/^ + TOo7±v^^ ( ±^Cif^'«(i) + T^u{ 2) + 2^±l'±^u(l) 


(4.17) 


The behavior of the energy E versus the angular momentum J in the case of co-rotating orbits is shown in Fig. 2 for 
selected values of the parameters. 

Finally, the 4-velocity U is given by Eq. (3.11), with normalization factor 


F = r± 


|l-h 7±i^±J^is-f 7 I 


V±V2 + 1 27± - 1 ) 



and angular velocity 


c 


c± 


1 + 


N 


C± y dcjicji 


{viS + I22S^) 


(4.18) 


(4.19) 


where 


= ^'u(i) , 1^2 = ^'u(2) ± 2 (1 - Cq) . (4.20) 

The relation between the timelike unit tangent vector U to the body center world line and the unit timelike vector u 
aligned with the 4-momentum thus reads 

U-u = y^l{l-CQ)MC,Kyity±. (4.21) 


the unit vector U± being already defined in Eq. (3.16). Hence, in general, U and P are not aligned unless Cq = 1, as 
discussed below. 


2. Weak field limit 


Let us study now the weak field limit of the above analysis. For convenience, we introduce the dimensionless 
quantities uq = M/tq and a = a/M. We may consider only the case of co-rotating orbits, the counter-rotating case 
simply following from the replacement a —^ —a. Every quantity is expanded up to a certain power of uq as follows 

X ^ Xq + aXa + sXg + o^Xad + ‘2-asXds + s^Xgg , (4.22) 

where terms of orders higher than the second in the background rotation parameter, as well as terms like as^ and 
a?s, are neglected. 

The weak field expansion of the conserved energy and angular momentum (4.17) are then found to be 


mo 

J 

moM 


+ 


—^0(2 -f 157/0)^^ 4- ”'^0(2 4- 23^0)0 


4-^0 
Mo 


/ 3 

27 2 135 3 

2835 4 

15309 gX 

0 5 

63 2 

405 3\ 

(1+2^0 + 

+ 

Mn -4 
128 ° 

256 “V 

— Olio 1 1 H" 2^0 H" 




-4 


81 .,3 4185 

“0 


27 2 81 o 

1 — 2mo -Un-Mn- 

“ 8 ° 8 ° 128 

3 3/2 A 9 159 2 1305 oN .. 

+ 2^0 (1 + 2 + “§““0 + 1 as 

23^ \ 

yEg 1 Mo 


Mq^^ ( 1 + 5mo 4- 


189 2 405 0 

—ao + —ao 


a 


+ 4^0 


C-Q- 




■0 j as 


ao + y(81 + 477CQ)M3 
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whereas the normalization factor (4.18) and the angular velocity (4.19) become 


r — 1 + -uq + ~ 16 ^°^ + s) 

+ tWq( 2 + 27uo)a^ + —itQ(2 + 19uo)as + —Uq[2Cq + 5(2 + Cq)uq\s^ + 0(uq) 
4 4 8 


M( = - ( a + -s j Mq + uy^a^ + -Mo'^^(1 + 2uo)as + 


C'q+(^-2Cq 


+ 0{ul), (4.24) 


respectively. In order to derive gauge invariant expressions, we express E and J in terms of the gauge-invariant 
dimensionless variable y = (MC)^/^, related to uq by 


Uo = 


^ ^ '®) (^1 - as - ^v^ICq - 2(1 -f CQ)y]s^ + 0{y^^y . (4.25) 


This yields 


E , 1 3 2 27 3 675 4 5/2 ^4^ 27 2\ . 5/2 f. 3 27 

— = + +16*' +128" -'■> (3++"+T" )“-" ('+ 2 "+- 


y^ I s 




CQ-h[l--CQ]y 


oiyy, 


J 1 / 2 /^ 3 27 2 135 3 2835 4 15309 /lO ^ 81 2 495 3 , 

= r+d 1 + j!/+ + ^ 1 /’+ ^!/‘+ » (y + 79 + - —»> I a 


21 


81 


1485 


^ 2 ^ 8 ^ 16^ 128 


y s + y 


26 335 


459 


+ y^'^^ (1 + —y + —y^ + -j-y" I a" 


9 


24 


3/2 4 25 2 81 3 \ _ 

-h ( 2 -y -f —y2 —y^ j as 


,,3/2 


Cg-4|l-lc«)!,-f (l-2co) „>-H(i_5c^)„3 




(4.26) 


3. Comparison with Ref. [36] 


Before investigating the general equatorial motion, let us show how the above analysis allows us to reproduce the 
results of Ref. [36]. As already stated, our quadrupole tensor reduces to the one adopted there to describe black holes 
ifwesetCQ = l. Since we have then 1 /= -1-0(s^) from Eq. (4.12) [see also Eq. (4.21)], this entails P oc C/-I-O(s^). 

In order to write our expressions for the energy and angular momentum in the same form as in Ref. [36] , we eliminate 
the dependence on the radial coordinate in favor of the angular velocity by inverting perturbatively Eq. (4.19): 


r = rc + —{-y^Mrc ± a)s -I- -^jAc - 2a(a=F \/ Mr^) - (1 - Cq) [Ac -I- 4a(a =F \/a^)]|s^ , 
where = y/M =F and Ac = A(r = Tc). Substituting next into Eq. (4.17) leads to 


E = 


mo 


ry^ — 2M^/r[ ± ay/M moM'^ 


3/4 r 


-k 


mg m3 


± 2ay/M 
1 


. 1 I /2 


T 


vMry =F a 


c9/4 


— 3My/r[ ± 2ay/M 


■1 1 / 2 ' 


9r-33/4 r 3/2 ^1 3/2 

rj - SA/^ ± 2ayjM 


^ 3/2 _ 4 - ^Qy/M 


- (1 - Cq) 


-C + T aaVM + liA + 3 ^ 3.) 

2v^ 2 rt' -3M^±2ay/M 


J = ± 


mo a/m Tc =F 2ay/Mr[ + < 


rngM r^(rc — 4M) -|- Ma^ ± ay/Mrd^rc — M) 


3/4 


± 2ay/M 


1/2 


9/4 

rc 


-md^±2ay/M 


1/2 


(4.27) 
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(a) 


(b) 


Figure 2. The behavior of the energy E versus the angular momentum J in the case of co-rotating circular orbits is shown in 
panel (a), for the fixed black-hole dimensionless spin parameter a/M = 0.5 and different values of s; in panel (b), for the fixed 
body spin s = 0.25 and different values of a/M\ Cq = 1 in both cases. Dashed curves in (b) correspond to the geodesic motion 
(s = 0). 


± 


mo 


M5/2 


2r. 


15/4 


-3M^±2aVM 


1/2 


3Ma^ ± 2aViWrI(3re - 2M) -h rl{2rc - 7M) 


- (1 - Cq) 


(2rc -h SM){rl - Mr^ + SM^) =f 4a/A^(3rc -h 2M) + lOr^a^ ± 


lly/Ma 

2Vr^ 


_^3y^6M3/2^(3M2 -h 4a2) =f a[a^{rc + 9M) -h 4M2(5rc -h 3M)] 
2 rl^ - ± 2a^/M 




(4.28) 


This exactly reproduces the results of Ref. [3G] when specialized to the case Cq = 1 [see their counterparts displayed 
in Eqs. (39a) and (39b) there, with in addition M = 1 = mo]. 


D. General equatorial motion 

1. Effective potential 


In order to discuss the general features of equatorial motion it is most useful to introduce appropriate effective 
potentials [48, 49]. The latter naturally arise when factoring the expression of (dr/dr)^ as a polynomial in the energy 
E of the test object. This factorisation takes the form 

= AV{E){E - R(+))(E - R(_)) + 0(s3), (4.29) 


dr Y 
dr / 


9 9 9 

7 COS a 


where the solutions for v and a in terms of the conserved energy and angular momentum are given by Eqs. (4.6)-(4.7). 
More precisely, we find 


A = 


- 2Ms6^{nY - ^fc(ue)(fc(Lie) - a{nr) 

+ (2 — CQ)(2Eff + Cgg )—2-b — CQ)Eng 3Effi 


/Vf2s2 

- Cq){E + N^J) 


{Eff + 2EgA{E + N<>’J) + 2 


N 

y/9H> 


-H,gJ 


(4.30) 
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and 


V^(±) = K) 


(±) 




(±) 




(±) 


(4.31) 


with 


yW _ 

Vq — 

= 

= 


T 

+ 


, N /- 

-N^J ± —— Jmlg^^ + J 2 , 
V5^ ^ 


MN 


M'^N f 


{a{ny + fc(Lie))J ± 0^{n) 


r ^^o9<P<P 4” 


± 


\/ mlg^^ + J2 ^ 

2^1 ± 2[0^(n)'-(a(n)’- + - CQH,g]J 

+ \J'mlg4,4, + .p [2fc(Lie)(fc(Lie) + 3a(n)’') - (4 - CQ)Eff 

W^mlN^/gy^ 


[fc(Lie)(fc(Lie) + ia^nf) - (3 - CQ)Eff - iEg^ 


ny66 “r J 


T^Q9(p(p 


(4 + CQ)Egg] 


(4.32) 


If E and J are kept fixed as s goes to zero, the leading order approximation of V{E) is obtained by setting exactly 
s = 0 in Eq. (4.30), which shows that 'P{E) is necessarily positive in the domain of validity of the small spin expansion. 
The solutions E = V(j-) are generalizations of the radial effective potentials for a test particle in a Kerr spacetime 
to the case of an extended body with spin-induced quadrupole moment. Their behavior as a function of the radial 
distance is shown in Fig. 3 for selected values of the parameters. The upper/lower branch corresponds to the -f/— 
sign in Eq. (4.32). 

On the other hand, since the equation (dr/dr)^ = 0 is quartic in E, in order to give a complete account of effective 
potentials we should also consider the solutions to the equation ViE) = 0, or equivalently. 


3aM^s^ 
2 ^ 


{2-Cq){E + N^J) 


a{E + N'^J) - 2^(r2 -h a^J 


= 0 . 


(4.33) 


If Cq < 2, no real solutions for the energy exist. If Cq = 2, then V^E) = 1, irrespective of E and J. If Cq > 2 (so 
excluding the black hole case Cq = I), the above equation admits two real solutions E = W(±), perturbatively in s, 
only if E scales as I/s, with 


1E(±) 


Too 

a 


± 


r n5/2 j 


+ 


s^/3{Cq- 2) yMJ moM 


0(s). 


(4.34) 


Note that the solutions kF(±) diverge in both limits a —>■ 0 and s —^ 0 for fixed values of the radial distance. 
Furthermore, for fixed values of the dimensionless spin parameters, they indefinitely grow for large r, exhibiting a 
monotonic behavior VF(±) ~ ±(r/M)^/^, so that there cannot exist circular orbits associated with them. In the 
following, we shall actually exclude the configurations for which E = 0(1/s) and ignore both potentials W(±y 


2. Circular orbits and IS CO 

Circular orbits correspond to the extremal points of Vy) and solving for J the resulting equation, = 0, provides 
the associated angular momentum. Now, as pointed out by Le Tiec et al. [50], the shift of the ISCO frequency due 
to the peturbation induced to the spacetime background metric by the particle itself is an important strong-held 
benchmark. Gravitational self-force theory has provided very accurate analytic predictions for it in the case of a 
spinless body in motion along a circular geodesic on a Schwarzschild background, at the hrst order in the symmetric 
mass ratio of the two objects. In the present situation, because of its spinning and quadrupolar structure, the body 
deviates from geodesic motion in the way described by the MPD model. As a result, the last stable circular orbit 
undergoes a shift in the frequency, made of terms proportional to the spin as well as the quadrupole. Measuring this 
effect can therefore provide relevant information on the structure of the body. Conversely, having information about 
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the spin and quadrupolar structure of the body allows one to make predictions on the frequency of circular motion 
and its deviation from the corresponding geodesic value. 

Stability requires that < 0. For a spinless object the latter condition boils down to 

- 6Mr - Sa^ ± 8a/Mr > 0 . (4.35) 


The equality can be analytically solved for the radius r^s of the marginally stable orbit (or ISCO) [51]: 


Kerr 
' ms 


M 


3 + ^2 T \/(3- Zi)(3 + Zi + 2 Z 2 ) 


(4.36) 


where the upper/lower sign refers to co/counter-rotating orbits, with 

= 1 + (1 - [(1 - + (1 + a)^/^j , Z 2 = + Zf . 


(4.37) 


The latter quantities are even functions of a satisfying Z 2 > Zi, Zi <?) and, for a = 0, = 3 = Z 2 . For small values 

of a, the Kerr ISCO radius may be expanded as 


„Kerr 
' nis 


M 


6t 


4^6, 


-a - a 

18 


13^6 

162 


a^- 


241 

1^ 


+ 0(a®). 


We recover the value = 6M in the Schwarzschild case (a = 0). 
For a spinning particle the ISCO is modified as follows: 


Kpt-t- ^ .n2 

^ISCO ^ms ~b '^^rrrs(l) 4“ ^ ^ms(2) ; 


(4.38) 


(4.39) 


where 


^ms(l) 

M 

^ms(2) 

M 


2^6 2^ , 17^6 ..2 

“ 9 “ 324 “ 

72 2 ^ 54 V 24 




Cq 


1 


451 

Is” 



± 


76 

1296 


4559 



7 + O(a^). 


(4.40) 


Of course, both the first and second order corrections to the ISCO can also be straightforwardly computed in the strong 
field regime. However, the corresponding expressions are quite long, so we prefer not to explicitly write them down. 
For instance, in the co-rotating case, for a = 0.5, we find fa 4.233, ryas{i)/M ~ —1.472 and r^s{ 2 )/^ ~ 0.194 

{Cq = 1), or r„is( 2 )/M k. 2.674 (Cq = 6). The resulting behavior of nsco as a function of the spin parameter s is 
shown in Fig. 4. 

The ISCO frequency of a spinning test object is then computed to be 


... , .-s/2 f. , 1 3 A 59 ,2 1 

MCisco = ±6 / |l±—(^-a+-sj+—a +-as + 


1 /97 


9 V 64 


-Cq]s 


leading to a fractional correction with respect to the spinless case 


^ISCO = 


Cisco 

/•Kerr 

MSCO 


-1 = ± 


76, 


■"+48“"+9 


1 /97 


64 


whereas the energy and angular momentum at the ISCO read 


^^isco 272 73, .. 272 

m^M 3 




1 


- O' H” —os - I — — C 0 ) s 

KA 01 I Q ^ / 


54 

11 


15 


216 V 8 


1 


—— d 4“ Qs F —T I f — xO/ 


27 


216 


24 


(4.41) 


(4.42) 


(4.43) 


The behavior of the ISCO frequency as a function of the spin parameter s is shown in Fig. 5 for two typical values 
of Cq, i.e., Cq = 1 (black hole) and Cq = 6 (neutron star). The corresponding curve is in general a parabola, which 
is concave up or down depending on whether the sign of the coefficient of is positive or negative. For instance, 
for the chosen values of the rotation parameter a = [0,0.5, 0.9], the change of concavity (from up to down) occurs 
at Cq ~ [1.516,1.487,1.455], respectively. Finally, Fig. 6 shows the behavior of the fractional correction to the 
ISCO frequency as a function of the spin parameter. Furthermore, restricting s to a given range of values yields the 
uncertainty associated with the ISCO position, angular velocity and shift. As an example, we list the three latter 
quantities below in Table I for selected values of the rotational parameter a = [0,0.1,0.3,0.5,0.7,0.9], spin parameter 
s = [—0.1,0, 0.1] and polarizability parameter Cq = [1,6], in the corotating case. 
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Figure 3. The behavior of the effective potential for radial motion is shown for the following choice of parameters: a/M = 0.5, 
Cq = 1, s = 0.25, and for different values of the body dimensionless angular momentum, J/{moM) = [1,2.5,4, 6]. The 
corresponding geodesic case (s = 0) is also shown for comparison (dashed curves). 



(a) 


(b) 


Figure 4. The radius of the ISCO as a function of the spin parameter is shown in the case of co-rotating orbits for different 
values of the black-hole dimensionless spin, a/M = [0, 0.5, 0.9], and for (a) Cq — 1, as well as (b) Cq = 6. 


3. Quasi-circular orbits 


Let us finally construct the quadratic-in-spin soiution to the MPD equations corresponding to a quasi-circuiar orbit, 
in the perturbative sense. The initiai conditions are chosen so that the worid iine of the extended body has the same 
starting point as the reference circuiar geodesic at radius r = tq for vanishing spin. We aiso require that the two 
worid iines are initiaiiy tangent. 

The orbit can be parametrized in a Kepierian-iike form as foiiows [52, 53] 


Stt 

— (t-to) = it - etsinit, 

r = Ur (1 — Cr cosir) ■. 

TT 


6 = 


2 ’ 


27r 

$ 


{4> — 4 >q) =2 arctan 



(4.44) 


where is some “semimajor axis”, et, Cr and are three different “eccentricities”, which wouid coincide in the 
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Figure 5. The behavior of the ISCO frequency as a function of the spin parameter is shown in the case of co-rotating orbits for 
different values of the black-hole spin, a/M = [0, 0.5, 0.9], and for (a) Cq = 1, as well as (b) Cq = 6. 



Figure 6. The behavior of the fractional correction to the ISCO frequency as a function of the spin parameter is shown in the 
case of co-rotating orbits for different values of the black-hole dimensionless spin, a/M = [0,0.5, 0.9], and for (a) Cq = 1, as 
well as (b) Cq — 6. 


Newtonian theory, while P and $ denote the periods of t and cj) motions, respectively (with an abuse of notation for 
P, not to be confused with the body’s 4-momentum). The quantities it-, ir and i^ are functions of the proper time 
parameter r on the orbit. They are conveniently expressed in terms of the dimensionless variable i = ri(ep)T, where 


^(ep) — l^±l 


6M 

^■o 


3a^ 

—^ ± aa/K 
^0 


1/2 


(4.45) 


denotes the well-known epicyclic frequency governing the radial perturbations of circular geodesics. The quantities 
it, ir, icj) and i play the role of eccentric anomalies. Note that, for geodesics, the above quantities reduce to Qr = tq, 
e-t = = Q, P = 27rr±, <i> = 2'nVL± and it = ir = i(p = i- For non-vanishing spin parameter s, the semimajor 

axis and the eccentricities for the orbit of the extended body turn out to be 


r — ^0 sRs S Cis , 


a, 


(4.46) 
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Table I. The numerical values for the modified ISCO position, angular velocity and shift are listed for selected values of s and 
a in the corotating case, for Cq = 1 (black hole), as well as Cq = 6 (neutron star). 






Cq 

= 1 






nsco/M 

Af Cisco 

5isco 

a \ 

-0.1 

0 

0.1 

-0.1 

0 

0.1 

-0.1 

0 

0.1 

0 

6.1643 

6 

5.8377 

0.0660 

0.0680 

0.0702 

-0.0300 

0 0.0312 

0.1 

5.8313 

5.6693 

5.5094 

0.0712 

0.0735 

0.0759 

-0.0311 

0 

0.0323 

0.3 

5.1351 

4.9786 

4.8248 

0.0847 0.0877 

0.0907 

-0.0335 

0 

0.0348 

0.5 

4.3824 4.2330 4.0876 

0.1046 

0.1086 

0.1127 

-0.0365 

0 

0.0376 

0.7 

3.5259 

3.3931 

3.3002 

0.1384 

0.1439 

0.1475 

-0.0379 

0 

0.0250 

0.9 

2.4018 

2.3209 

2.2427 

0.2194 

0.2254 

0.2317 

-0.0267 

0 

0.0277 





Cq- 

= 6 






nsco/M 

MCisco 

disco 

\ ^ 
a \ 

-0.1 

0 

0.1 

-0.1 

0 

0.1 

-0.1 

0 

0.1 

0 

6.1893 

6 

5.8627 

0.0656 

0.0680 

0.0698 

-0.0356 

0 

0.0256 

0.1 

5.8565 

5.6693 

5.5346 

0.0708 

0.0735 

0.0755 

-0.0370 

0 

0.0264 

0.3 

5.1604 

4.9786 

4.8501 

0.0841 

0.0877 

0.0901 

-0.0402 

0 

0.0281 

0.5 

4.4065 

4.2330 4.1122 

0.1038 

0.1086 

0.1119 

-0.0438 

0 

0.0301 

0.7 

3.5438 

3.3931 

3.2815 

0.1375 

0.1439 

0.1486 

-0.0446 

0 

0.0330 

0.9 

2.4094 

2.3209 

2.2504 

0.2184 0.2254 

0.2306 

-0.0313 

0 

0.0230 
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respectively. As for the periods of t and </> motions, they may be written as 


P = 2-k 


r± 




(ep) 


1-1-0 -\i(r)* I Djig ,2 

1 =F 2'y±v±— -V) 's + ——s 


PI 


(ep) 


r± 


$ = 27r- 
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(ep) . 


1-1-0 1 jP) ~ I -2 

lT2j±iy±- -—Vi 's+—s 

“(ep) < 5 ± ll± 


There remains to display the three eccentric anomaly parameters, £t, ir and as functions of l\ 

(■t = £ + — [Il 2 sf^(ep) sin 2 t + Dsgi cos £], 


£r = £ + 


£th — £ 


^sf^(ep) 


[2C'2sf^(ep) sinf-Csj^], 




^(ep) ( I sin 2£ -b E 3 g£cos£ 


(4.47) 


(4.48) 


(4.49) 


The various coefficients entering the above formulas are listed in Appendix E. Notice that the parametrization (4.44) 
of the orbit is different from the quasi-Keplerian one, used in Refs. [52, 53], due to the presence of different parameters 
representing the t, r and </> motions instead of a single eccentric anomaly. The two forms agree to first order in spin 
only, as shown in Ref. [33]. 

At this stage, we can derive explicit expressions for the conserved energy and angular momentum (2.15): 

^ ^ ^ ^ ^ J ^ 

E = — = E± + sEg + s'^Egg , J = — = L± -b sJg + s'^Jgg , 

mo mo 


(4.50) 
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with 


Es = Mn± 


M 

ro 


T OjCk 


Eg C±'^S - T 


MCk 

r± ’ 


(4.51) 


and 


hs = -V^llvP 


(1-Cq)M%p) 


— 


7±i/± 


7712 r 

mo 


Ess - C±Jgg = 2N-t±{V^pf + , 

7770 r± 


(4.52) 


where E± and L± are the energy and the azimuthal angular momentum per unit (bare) mass for a circular geodesic, 
as given by Eqs. (3.18), the parameter m 2 representing the spin correction (4.8) to the mass of the body [see also 
Eq. (4.13)]. 

In the weak held limit the previous expressions become 


—2E = — 2{E — 1) = mq ( 1 ~ 7^0 —Uq —— 


27 


675 


64 


9 5/2 / 1 , 9 , 135 2 1 - 

2 uo ( 1 + 2^0 + “'“0 ) a 


— 2uq^^ ( 1 + + TT^o 1 s — Wq ( 1 + ^uo ) a‘‘ + 2 uq ( 1 + 1 as 


27 


15 


5 

12 


Cq +6(1 + —Cq ] Uq + 0{u, 


1). 


J=— - 

"“0 


. 3 27 2 

1 E —Uq H- ^Uq + 


, . 1 3 2 27 . 

4 “ 1 — nUQ + —Uq + 777U0 


16 


135 

675 

' m 


2835 

128 


63 


Uq — Suq 1 + -Uq + —Uq + ——Uq a 


405 


16 


Mq I S + u] 


3/2 


0 


- 21 + 


5/2 

0 


14 


^ + 42*^^ + T7 ( + O^Q ) ^0 + 


15 


1 + 5mo + 
3 


189 


405 


Uq + -^Wn I — Sul.^'^aS 


111+2^0 


s^ + 0(mo) ) 


(4.53) 


for the co-rotating case. We list below the orbital elements (4.46)-(4.48) expressed in terms of the gauge invariant 
quantities E and J, related by 


J = (-2E)-i/2 \ 1 + ^(_2E) + ^(-2E)^ + §§(-2^)' - 2(-2E)3/2 


l + ^(-2E) + ^i-2Ef 


+ (_2+)1/2 


1_ 3(_2g)_ 9(_2^)2_|L(_2^)3 


1 


+ i-2Ef 


l + ^{-2E) 


as + ^{-2Ef 


s+-{-2Ey 


7 


69 ' 

l + ^(- 2 E) 


Cq - 3 15 - -Cq (-2C) 


52 + 0[{-2Ef] , (4.54) 


so that we may write them with the help of a single parameter, e.g., the energy parameter (—2E) = X. We hnd 


M 


= X 


-1 


1 _ ^ 
4 16 64 256 




63 


-X [1 + —X + 183X2 ]a^ -X [1 + —X + 213X 


79 


as 




Cq - ( 36 + 7 CQ 1 X - 3 ( 86 - icQ 1 X 


+0(X^/2), 


et 


= -6X5/2 / (^1 + s - Xi /2 (1 + —X^ as - ixi /2 
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Cq + ( 1 + ^Cq ) X 


s^+OiX^)\, 


= 3 X 5/2 I s - xC 2 [1 + ^x^ as - Ix^/^ 


= 6X5/2 ) C ^ \ - _ 


X5/2 (^1 + ^X^ as - ix5/2 


Cq- 10--CQ X 


Cq-(4--Cq)X 


+ 0(X2) , 


s^ + 0{X^)} . 


(4.55) 
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At last, the periods of t and (j) motions are given by 


— — = X-3/2 11 + —X + —X^ + - 12X5/2 f ^ ^ - _ 3^3/2^^ + 6X + 39X2)s 

M2tt \ 8 128 1024 \ A J v t t ; 


+llX5a2 + 3X2(1 + ^X^ as + ^X2 




1 + ?Cq)X 


s2 + 0(X4) , 


$ 63 , 405 . / 93 \ o/n / 57 \ 

— = 1 + 3X + —X2 + —X5 - 4X^/2 1 + —X a - GX^/^ 1 + _X 

27r 44 \ 8 J \ ^ J 


73 


+ -X^ 1 + —X + 6X^ 1 + —X as + 3X 


39 


Cq + -(13 + 11Cq)X 


s2 + 0(X^/2) . 


(4.56) 


V. CONCLUDING REMARKS 


We have investigated finite-size effects on the motion of extended test bodies, in the equatorial plane of a Kerr 
spacetime, within the framework of the Mathisson-Papapetrou-Dixon model up to the quadrupolar order. In general, 
the quadrupole tensor shares the same symmetries as the Riemann tensor and is completely specified by two symmetric, 
trace-free spatial tensors, i.e., the mass quadrupole (electric) and the current quadrupole (magnetic) tensors, whose 
role has been investigated in previous works [33, 34]. Here we have considered the rotational deformation induced by 
a quadrupole tensor of the electric-type only, taken to be proportional to the trace-free part of the square of the spin 
tensor, with a constant proportionality parameter which may be regarded as the polarizability of the object. This 
allows us to treat on an equal footing the cases of black holes and neutron stars, so generalizing previous works. 

The general features of equatorial motion have been discussed through the analysis of the associated radial effective 
potentials. We have obtained their generalization from the well-known case of a co/counter-rotating test monopole 
particle in a Kerr spacetime to that of an extended test body with spin-induced quadrupole moment. We have also 
evaluated the correction to the ISCO due to spin and the corresponding frequency, which is an important observable 
in gravitational-wave astronomy. The presence of spin corrections introduce an uncertainty on the values of the 
corresponding quantities for structureless particles. On the other hand, those features can be used to determine 
whether the small object is endowed with a spin, by performing an adequate parameter estimation in the context of 
gravitational-wave detection. 

The dynamics of the system have been studied not only qualitatively, but also quantitatively. In fact, neglecting 
terms in the MPD equations that are of third order in spin or higher allowed us to solve the problem in a full 
analytic way. Initial conditions have been chosen so that the tangent vector to the orbit of the extended body be 
initially tangent to the 4-velocity of a timelike spatially circular geodesic, taken as the reference trajectory. We have 
obtained the “perturbative” solution to second order in spin in the following two cases: (i) when the trajectory of the 
extended body remains circular with spin-dependent frequency, (ii) when it deviates from circular motion because of 
the combined effects of both the spin-curvature and quadrupole-curvature couplings (i.e., when the orbit is “quasi¬ 
circular”). The tangent vector to the orbit and the unit timelike vector aligned with the 4-momentum are in general 
distinct. However, there exists a special value of the polarizability constant, which corresponds to the black hole case, 
such that they are aligned not only initially, but all along the (circular) trajectory of the extended body. This is no 
longer true for neutron stars, an interesting fact which seems to have never been pointed out before. For quasi-circular 
orbits, we have explicitly written down the solution in a Keplerian-like form, by introducing the temporal, radial and 
azimuthal eccentricities of the orbit, as well as the associated periods and frequencies. We have also computed the 
spin-induced shift of the conserved energy and angular momentum, in a gauge-invariant way. All orbital elements 
have been expanded in the weak field and slow motion limit, in a more suitable form to be compared with the existing 
post-Newtonian literature. 
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Appendix A: 1 + 3 decomposition of the quadrupole tensor 

Let us consider an orthonormal frame adapted to an observer family characterized by the 4-velocity u (with nor¬ 
malization u ■ u = —1), say eo = u and {cq}, a = 1,2,3, a triad of three unit spatial vectors orthogonal to u. 
We shall introduce the compact notation Xap..u°‘ = for tensor contraction. In addition, we shall denote by 

*Xpa... = paXap... and X* ^ the left and right dual of a tensor, respectively. The standard 1-1-3 

decomposition of the Riemann tensor in terms of its electric (spatial and symmetric) part E(u), its magnetic (spatial 
and tracefree) part H(u), and its mixed (spatial and symmetric) part F{u), defined by 

E{u)afi = Raufiu , H[u)al3 = ~RauPu i F{u)ap = J (-^1) 

respectively, leads to the identification of the 20 original independent components: 6 in E{u), 8 in H(u) and 6 in 

Fiu)- ^ ^ 

Similarly, since the algebraic symmetries of the quadrupole JafijS are the same as for Rap-yS, one can decompose 
the former quantity in terms of the associated tensors 

Q{u)aP = JauPu , LV(M)ct/3 = ~JauPu J Xi{u)ap = J^uPu ■ (-^2) 

In so doing we identify its electric (spatial and symmetric) part Q{u), with 6 independent components, its magnetic 
(spatial and tracefree) part W(ti), with 8 independent components, and its mixed (spatial and symmetric) part A4(u), 
with 6 independent components. However, J enters the MPD dynamics only in certain combinations, through specific 
contractions with the Riemann tensor or its derivative. Hence, the number of effective components needed is reduced 
by half, as shown in detail in Refs. [32, 33]. The proof requires the replacement of the mixed part A4{u) by a new 
tensor X (u) (with the same symmetries), according to 

M{u) = Q{u) + X{u), (A3) 

as well as the decomposition of both X{u) and yV{u) in terms of their STF and pure-trace parts, 

X{u) = Xiuf'^^ + i[TrA’(u)]P(u), W(m) = W(m)®^^ -f \[TryV{u)]P{u ), (A4) 

o o 

where [H(m)]“^ = <5“^ -I- u^'up denotes the projector to the hyperplane orthogonal to u. Inserting the resulting 
expression for J into the equations of motion then cancels the contribution of Q(u), which yields the following 
“effective” representation of the quadrupole tensor (valid only in the context of the MPD model): 

^“^75 = + 2m[“[W(u)]®'^^^1<^77(u)'^.^ 5 -h 2uY^[W{u)f'^^s-^a-niuYap , (AS) 

with r]{u)ap~f = u^ripap-y dehning the space 3-volume form (see section H). Summarizing, in basis components, we can 
write 

= [A(u)STF]-, , 

J°\c = = [W(u)®t^]*(“)%c, 

r\d = = r(“)[A(u)STF]A-)]-b^, . (A6) 

For convenience, we actually use the notation: 

X{uf'^^ = X{u ), W(it)®'^^ = W(u). (A7) 


Appendix B: 1 + 3 decomposition of the MPD equations 

It is useful to perform a 1 -|- 3 splitting with respect to U of the MPD equations (2.1)-(2.2). A key observation is 
that the force term on the right-hand side of the first equation (2.1) is not spatial for the comoving observer with 
4-velocity U, since F'(spin) -11 = 0 whereas F’(quad) ■ U 0. Recalling that the operator P{U)°‘^ = S°‘p + U°‘Up 
represents a projector perpendicularly to U, we see that 

= ^(':pin) + [^^(C^)^(quad)]^ + \u^ 


dr 
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^ jajB'yS 

dr 


(Bl) 


where the force 

^(t^)(tot) ^ ^(Vn) + [^(C^)rquad)]'^ = [^(C/) «pin) + J"(quad))]'^ (B2) 

is orthogonal to U and the mass correction mj has been defined in Eq. (2.13). In a second stage, we get from Eq. (Bl) 

or, equivalently. 


(B3) 


DP^ dmj 1 D 0 g 
Un— ^ ; = — J • 


dr 


dr 


PiU) 


DPM 

dr 


6 dr 


= nu)^tot) 


(B4) 


which follows from projecting Eq. (B3) along U and perpendicularly to U. 


Appendix C: ZAMO relevant quantities 


We list below the non-vanishing components of the electric and magnetic parts of the Riemann tensor, as well as 
the relevant kinematical quantities as measured by ZAMOs and evaluated in the equatorial plane. 

The radial components of the acceleration and expansion vectors are given by 


i{nY = 


M {r^ + ~ ‘icL^Mr 


{nf=-- 


aM{3r'^ + a^) 


r'^^/A + a^r + 20^]^ ’ {r^ + a'^r + 2a'^M) 

respectively. The expressions for the radial components of the curvature vector are 

Mr - 


(Cl) 


At(r, nY = 






(r^ — a?M)yYA 
r^(r3-I-a^r-f 2a^M) 


(C2) 


Finally, the nontrivial components of the electric and magnetic parts of the Riemann tensor with respect to ZAMOs 
read 

M{2r'^+ ^"^ 0 ?— 2a?Mr+ 'iaY ^ j? t? jr ^ 

’ ^§§ ~ ~ ’ 


r^(r3 _ l _ q2j, _|_ 2a^M) 

3Ma{r'^ + a'^)'/A 
r9 r'^{r^ + a^r + 20^ M) 

In the limit of vanishing rotation parameter (a —?► 0), the previous quantities reduce to 

M 


a{nY = 


Nr^ ’ 


/ ^ , N ^ 2M 

(?l) — 0 , ^(Lie) — ; Ef-f — — 




(C3) 


(C4) 


with N = y/i — 2M /r. 


Appendix D: Frame components of both spin and quadrupole terms 

We list below the explicit expressions for the components of both spin and quadrupole terms with respect to the 
frame adapted to u. 

The spin force is given by Eq. (3.25) with 

^(spin) = silu {vu cos 2auEf.f + [v cos(a„ - a) + Vu cos^ au]Egg -f sin [1 -|- Wu cos(au - a)]H.g } , 
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^(spin) = -sgn(z/„)s 77 ^ {z/„[sin2au - vvu sin(Q!u + a)]Ef.f + [i/sin(Q;u - a) + Vu cosa„(sinau - Wu sina)]i?gg 

-[cosau(l + Wu cos(au - a)) - 2wuC0sa)]H~g] , (Dl) 

the remaining component following from the condition -/^(spin) -17 = 0, i.e., 

7 u(l - Wu cos(a„ - a))-F(°pin) = sm{au - a) + sgii{uu)F^^p-^^'yu{-i^u + i^cos(a„ - a)). (D2) 

On the other hand, the spin quantity I?(spin) takes the form of Eq. (3.26), with 

f(spin)(u) = mj [j^sin(a„ - a)uj'^ + sgn(^'„) 7 „(^ cos(a„ - a) - Vu)uJ^] ■ (D3) 

Concerning the quadrupole contributions, one gets for the components of the force (3.25): 


rpl 

^ (quad) 


/7’2 

^ (quad) 


1 ^2 r 

= ~^~^Q \ [biVu sin3a„ + 2(64 cos 2 au + 65 )] - 

1 52 r 

= -sgn{vu)-.—CQ^u{ bi'-fli'lcosSau + 2 [(oi - b4)^l - ai] r'„sin2au 
4 m 


(61 - 262 ) 7 ^ - bi + - {b2 + bs) 


sm a, 


(&i + 203 ) 7 ^ -bi- -{ci + C 2 ) 


cos a. 


with 


ai = {Fgg + 2Fff)0^{ny + H.ga{ny , 
a2 = (2£;,-^ + Frr)a{ny - H^gd^inf , 

bi = —2{Fgg + 2Ff.f)k{i^{e) + (2E’gg + Ff.f)a(ny — AFl.g9^{ny + —dfFg 


99 ’ 


62 = -{Fgg + 2£’rf)fc(Lie) + ‘^H^g0^{ny - -dfFgg , 

63 = -2{Fgg + 2Ffr)k(^ue) - ‘^H-g9^{ny , 

64 = H.gk(^i^ig'f — 2>Fgg9^(n) + dfH.g , 

b5 = i?fefc(Lie) - (Fgg + 2Fff)9^{ny - dfH.g , 


(D4) 


(D5) 


as well as ci = C 3 + 2 a 2 , C 2 = C 3 — 62 and C 3 = [762 + 4 a 2 + 6(61 — 263 )]/ 5 , whereas E|Pjuad) obtained by 

requiring that the coordinate component F^^^uad) t = 0 vanishes, which implies 


0 = 7«(^’(°quad) + sgn(z/„)z/„F(2qu,,d)) + 




N 


^(quad) COSa„ - 7„ sinau(l^uE(° ,^ad) + Sgn('^«)^(quad)) 


(D 6 ) 


Finally, the torque tensor E)[(^^ad)’ whose structure is displayed in Eqs. (3.26)-(3.27), is such that 

gS 

^(quad)i ~ COS \yu sinG;.^ 4 “ ^§0) -^f^] ’ 

s 2 

'^^(quad )2 = sgn(r'a)— 0 ( 37 ^ [uu cos2auFff + ^'„(1 + cos^ au)Fgg + sina„(l + vi)H.g] . 


(D7) 


1. The Schwarzschild limit 

We list below the corresponding expressions of both spin and quadrupole terms in the limit of vanishing Kerr spin 
parameter. 

The spin force (3.25) becomes 
M f 

-^(spin) = 37„:/„sina„(^cosa - i/„cosaa)u 

+ t[j/^(l - 3 cos2q!u) + 2 i/cos(Ga - a)]ei 
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, 3MNs^ 2 . 

^(quad) =- 


1 - + 5cos2a„) 


3MN 


^(quad) = -sgn(i/^)^^—CQ7„cosa„ 
Finally, for the torque term as shown in Eq. (3.26), we have 


1 + - 5cos2q!„) 


3M s 


'(quad) = --^—CQ-fuVu sin 2 a„[w^ - sgn(^'„) 7 „ tana„a;^]. 
2r'^ m 


- -sgn(z/„) 7„[(2 + v‘l)vs\n{au - a) + ^Vu{vvu sin(Q;u + a) - sin 2 a„)]e 2 , 

whereas the spin quantity I?(spin) is still given by Eq. (3.26) with components (D3). 

The quadrupole force (3.25) reduces to 

i^(quad) = -F(Vad)ei + -^'(quad) [-Sgn(l^u)j^«U + 63 ] , 

where 7 „[—sgn(i^u)t'iiM + 63 ] represents a unitary and spacelike vector orthogonal to n and 

1 


(D 8 ) 


(D9) 


(DIO) 


(Dll) 


Appendix E: Quasi-circular orbits: coefficients 

We list below the various coefficients entering the quasi-circular orbit solution (4.44): 

,ll^± Ck ,,(r) 
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TV 02 a 

(ep) 
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(r) o 7±C± ^ a j 

= -^—Fi - K7T- — T ^oCk , 

ro^(ep) \ro 


(El) 


(E2) 


and 


Cis = 7± 


C2s = 7± 


^ 0 ^(ep) 

2 ro^^(ep) 


|i^±(Bi^D(ep) - Sas) +7±(vf^)2 {K{r,ny + fc(Lie))} , 


I^±B2s - 


7± 


2 Q 


(ep) 


(K(r,n)'' + /c(Lie)) (vf^)2 


^ ^ D 

Das = -7±i^±—- Bss , 

T’Oiqep) 


)(r) 


- 2 B„) + 

(ep) (ep) 


+TWD(ep)7ivi’') 


1 T ^2 ((3^± “ 4)A:(Lie) ± 6i/±Cif) 


27±^±(vr)^ 


(ep) 


7~\ _ __ 2 2 j-y 

3s - 


(ep) 

7i 77,(A 72 [ 7±4 „{r) , 1 _ ^±Ck r,.,„3 


4TVQ2 

(ep) 


(vD^ 


I 4TWQ(ep) 


V^ + lT 


Q 2 

(ep) 


[(3^1 - 4)A:(Lie) ± 6 j7±Ck] 


I 




Da, = t27^477^^3s . 


■ A 02 
^''“(ep) 


D4, = —Q(ep)(Dis + 2D3,) — Das, 


Dis — C±17i, — 


7± 


^ih y^^00fl(ep) 


(vfV, 


-£'2.5 — C±-£l 2 s 


7± 






(D ',2 



























25 


Ess — C±D3 s ■. 


Eis — C±Eis + 
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